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NASA TECHNICAL MEMORANDUM 

HIGHLY AUTOMATED OPTICAL CHARACTERIZATION WITH FTIR SPECTROMETRY 


INTRODUCTION 

When characterizing crystals for infrared sensor applications, it is im- 
portant to establish quantitatively the spatial uniformity of properties 
affecting the response to the infrared phonons. Optical characterization 
using infrared transmission edge mapping is of particular interest for alloy 
semiconductors such as Hgj_ x Cd x Te because the transmission edge is directly 
related to the composition of the alloy. While the method of obtaining opti- 
cal spatial characteristics by using a small aperture in the optical path and 
moving the specimen behind the aperture has been employed previously (see 
reference [1]), a facility with advanced capabilities has been developed. The 
capabilities of the facility in addition to its associated hardware and 
software will be discussed. Moreover, the facility will be shown to be useful 
in assessing the correlation between the crystal growth parameters and the 
compositional redistribution during the solidification of narrow band-gap 
semiconductors. [2] In this laboratory, I I -VI alloy semiconductors are grown 
from the melt by directional solidification [2-4] and quenching. [5,6] Several 
results of the transmission edge mapping technique have been helpful in 
clarifying the fluid flow phenomena occurring during these crystal growth pro- 
cedures, and the internal temperature gradients [7] and effective diffusion 
coefficients [8-11] have been estimated in the alloy melts. Some of these 
results will be demonstrated with spatial maps of optical properties such as 
composition. 


HARDWARE 

The equipment for the facility is a Mattson Instruments, Inc. Sirius 100 
Fourier Transform Infrared (FTIR) spectrometer system. The spectrometer in- 
cludes an x-y stage which allows for desired areas of the sample to be scanned 
with minimal operator input. The sample is placed in front of a pinhole aper- 
ture, and the optical beam is focussed to an area 250 ym in diameter at the 
aperture to maximize the power to the sample. The metal foil containing the 
aperture does not actually touch the sample, but stays a few micrometers ahead 
of the sample in the optical path. [12] 

The spectrometer is currently set up to operate in the mid- infrared 
region (2-15 ym) to accumulate optical transmission data but may be set up for 
other capabilities. In general, the system may be used to acquire reflection 
as well as transmission data through near-, mid-, and far-infrared regions. 
Useful spectra have been obtained in the mid- infrared range with apertures as 
small as 25 ym in diameter; however, the results shown in this report were all 
obtained with a 100 ym aperture which gave the optimum compromise between spa- 
tial resolution and the time required to obtain spectra with satisfactory 
signal -to-noise ratios. The beam of the signal sent through the sample is 
focussed and detected by a bulls-eye mercury cadmium telluride detector which 
is cooled at liquid nitrogen temperature to detect photon excitations. Since 
a small aperture is used, all signals fall on the central part of the detec- 
tor. Only the central part is used so that signal -to-noise is increased over 
what would be obtained with the entire detector. [12] 


The detector, beam splitter, and infrared source may be changed as 
necessary to match the spectral range of the spectrometer to the infrared 
response of the material. For the samples in this report, the infrared source 
was a globar which provides near-blackbody radiation, and the beam splitter 
was potassium bromide (KBr). The bench is equipped with a movable stage which 
can be moved in steps in the vertical plane by electronic stepper motors lo- 
cated at the base of the stage. One step is 50 ym, and it is the smallest 
increment the stage may move. The sample on the stage may be positioned 
throughout an area of 2.54 x 2.54 cm . 


SOFTWARE 

There are several software capabilities provided by the Mattson manufac- 
turer which allow for generating a spectrum and then manipulating the spectrum 
to obtain specific characteristics, such as finding peaks in a spectrum, or 
smoothing, integrating, or differentiating a spectrum. The system makes use 
of a header file which includes many parameter settings for proper scanning 
and spectral collection. Some of the parameters include number of scans taken 
per point (nscans), sample shape (shape), sample size (radius for a circle or 
x-y coordinates for rectangular wafers), detector type (det), etc. (see sample 
header listing, Table 1). 

Much of the software has been generated by the authors, in-house, to 
supplement and enhance the capabilities of the system. The spectrometers and 
terminals are run through a main computer by the UNIX (a trademark of AT&T 
Bell Laboratories) operating system. While most of the author-generated 
software was written in the C programming language, the UNIX system generally 
makes it convenient to handle unusual conditions. For example, "noisy" 
spectra can be smoothed and re-analyzed. Or, if the spectra contain impurity 
peaks such as carbon dioxide near 2350 cm -1 (4.26 y m) because of inadequate 
nitrogen purging, these peaks can easily be replaced by a straight line 
segment connecting any two input wavenumber points (using software provided by 
the spectrometer manufacturer) before the spectra are analyzed. [12] The 
software generated by the authors uses additions to the header files used by 
the standard spectrometer software (see Table 1). Software has been created 
to read the header file additions and move the sample with the x-y stage 
according to specifications. Once a spectrum is produced by performing a 
Fourier transformation, the cut-on wavelength is found by the following 
method . 

Parameters are added to the header file which include transpl, transp2, 
opaquel, and opaque2. The first pair of parameters is set in one of two ways, 
described below, and is used to determine a transmittance value characteristic 
of each spectrum which is called transp. The second pair is set in a range 
where there is no transmittance (see 0i and Op in Figure 1) and is used to 
determine a parameter called opaque, this method of designating opaque 
regions is used instead of merely assuming zero transmittance in the opaque 
spectral region to account for spectra with shifted basel ines. [12] The algo- 
rithm locates the portion of a spectrum with transmittance greater than 25% 
and less than 75% of the span between opaque and transp and then fits a 
straight line to the points in that region using the method of least squares. 
The intercept of this line with the wavenumber axis is taken as the cut-on 
wavelength. 
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Transp is dependent on the two parameters, transpl and transp2, which are 
set in a range where the spectra exhibit transmission (see Ti and Figure 
1). The choice of transpl and transp2 is determined by two different methods 
depending on the amount of free carrier absorption present in the sample. The 
first case involves samples with low free carrier absorption for which the 
transmittance is essentially constant over a substantial range beyond the 
cut-on. The range designated by transpl and transp2 is simply averaged and this 
value is used as the value of transp. In the case of high free carrier absorp- 
tion, in which the transmittance rises to a peak and then falls off im- 
mediately at lower wavenumbers, the wavenumber corresponding to the maximum 
transmission between transpl and transp2 is located (see T max in Figure 1). 

Then, beginning at the wavenumber for T (nax and extending to Tower wavenumbers 
over a range set by another additional header parameter, dT, the transmittance 
is averaged to obtain the value of transp (see the cross-hatched region in 
Figure 1). 

The sample is mapped in a manner which depends on its shape. For a rec- 
tangular sample, the mapping is started in one corner and proceeds back and 
forth horizontally with a vertical step at the completion of each horizontal 
line. For circular samples, mapping begins in the center and one half is 
mapped. After this half is completed, the sample is returned to the center 
and the second half is mapped. Mapping may be specified to a small region, 
such as a single point or a line, or other isolated area. For a sample with 
at least 15% transmittance, the time required for each point is about one 
minute including interferogram acquisition and the fast Fourier transform 
(FFT). Longer times are required for samples with less transmittance. The 
FFT is computed, at least in part, while the next interferogram is accumu- 
lated. To minimize data storage requirements, only the processed spectra are 
saved. [12] 

The resolution normally used to collect spectra is 8 or 16 cm"*. 

Spectral data points computed by the FFT are equidistant in energy, and each 
value of resolution above corresponds to a data point spacing of 3.86 or 7.72 
cm' 1 , respectively. [12] 

After samples are mapped and analyzed for cut-on wavelength, the results 
may be used to reveal the characteristics of the sample. For Hgj. x Cd x Te 
samples there is a known correlation between cut-on wavelength and the com- 
position, x, of the sample. Several spatial composition maps have been gener- 
ated for this alloy from the transmission spectra obtained. For materials 
such as Hgj_ x Zn x Se, the correlation between cut-on wavelength and composition 
is not yet as well established. Instead, spatial plots of wavenumber charac- 
teristics have been generated to show the spatial trends in wavenumber (or 
energy gap) of the sample. Several characteristics of a spectrum are stored 
in a file which includes the spectrum-file name (name identifying the sample 
and the coordinates at which the spectrum was taken), cut-on wavenumber value, 
percent transmission of the spectrum, and a composition value (if a correla- 
tion exists between wavelength and composition). This file is then used as 
the input for programs used to create spatial maps. Other plots may be gener- 
ated to show more clearly the characteristics of a sample caused by growth 
conditions such as the solid-liquid interface shape during growth and thermal 
asymmetries in the growth furnace. 

Only spectra with a transp value above a particular minimum (set by the 
additional header parameter transmin) are catalogued into the file used to 
create the maps. The minimum is usually set to eliminate spurious contribu- 
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tions from opaque areas of the sample or parts of the sample holder that were 
inadvertently included. Some low but detectable transmission may be rescanned 
using more scans to improve the signal -to-noise ratio and thus the spectra may 
be easier to analyze. Occasionally, in samples with large compositional 
variations (large variations in the cut-on wavenumber) and with high free 
carrier absorption, the values of Tj and T 2 must be reset for different parts 
of the sample. This problem may be alleviated by choosing wider limits on Tj 
and To at the expense of slightly longer analysis time. Also, these samples 
usually show very low transmittance near the edges which means the associated 
spectra will have low signal -to-noise ratios. This particular problem has 
been addressed by attempting to check the signal -to-noise ratio and then, as 
stated above, taking additional scans if the ratio is below a preselected 
level. When very inhomogeneous samples are to be characterized, these 
problems require additional efforts for proper analysis which include manual 
determination of cut-on wavelengths and spectral manipulation. Computer 
programs will be modified to accommodate extreme cases. 


RESULTS 

The Mattson spectrometer has been used to characterize many bulk samples 
of I I -VI semiconducting alloys including Hgi_ x Cd x Te, Hgj_ x Cd x Se, Hgj. x Zn x Te, 
and Hgi x Zn x Se grown in-house either by unidirectional solidification or cast- 
ing. Several epitaxially grown and cast grown Hgj. x Cd x Te samples from other 
labs have also been analyzed. 

Figure 2 shows a spatial map of CdTe composition, x, of a 5 mm diameter 
cross section of a Hgj. x Cd x Te ingot grown by unidirectional solidification. 

The data were taken on a grid 4.4 mm in diameter. The center is shown to be 
Hg-rich which is typical of a directionally solidified Hg 1 . x Cd x Te sample. In 
addition, the map displays radial asymmetry which is caused by asymmetries in 
the thermal field. The Hg-rich center indicates an interface shape concave 
toward the solid. These illustrations provided by the spectrometer system es- 
tablish the correlation between the compositional distribution of the grown 
crystal and the growth parameters and will provide useful information for the 
control of the solid-liquid interface shape and radial thermal asymmetry. 

Figure 3 is a spatial map of an axially cut wafer of Hgj_ x Cd x Se. The 
slab dimensions are 5 x 19.3 mm and the data were taken over an area of 4.5 x 
19 mm. During growth, the furnace translation was stopped for 31 hours. [7] 

The map shows buildup of Cd and the flattening of the isoconcentration 
surfaces as the actual growth rate decreased; also, the Cd-rich center typical 
of unidirectionally solidified Hgj. x Cd x Se alloys is apparent. From results 
such as this, the mechanics of growth processes may be better understood and 
controlled to get desired results. 

As mentioned previously, a direct correlation between the transmission 
edge cut-on wavelength and composition has not been established for all 
materials. Figure 4 shows a spatial map of wavenumber for a 5 mm cross 
section of a unidirectionally solidified Hgj_ x Zn x Se ingot. The data were 
taken on a 4.4 mm diameter grid. Transmittance values may also be displayed 
by spatial maps for all samples if desired. 

Other figures may be generated including three-dimensional plots and 
contour maps of alloys. These figures may be generated by DISSPLA (a product 
of Integrated Software Systems Corporation) which is available on the Marshall 
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Space Flight Center Engineering Analysis and Data System (EADS). Mapping data 
can be uploaded onto EADS. Figure 5 shows a three-dimensional plot of a 
Hgj. x Cd x Te cross section (same cross section as in Figure 2). The plot shows 
the relation between the composition and the location on the sample which is 
related to the interface shape of the material. 

A contour map of a directionally solidified Hgj_ x Cd x Se sample is shown 
in Figure 6. This map shows the Cd-rich middle typical of these samples and 
illustrates the trends in composition and uniformity radially on the sample 
cross section. Note, however, that the composition at the edges of the con- 
tour map are not necessarily the true compositions in those regions unless the 
sample happened to have those compositions at those points. For both the 
three-dimensional plots and the contour maps, the average composition is set 
at the edges. 

Combining several average compositions from different cross sections 
along a directionally solidified ingot, the growth effects may be assessed 
axially. Figure 7 displays the change in CdTe composition, x, along the 
growth axis of a 12.5 cm Hgj. x Cd x Te ingot which was unidirectionally 
solidified using a transverse magnetic field. [13] Note how well the FTIR data 
compare with the data taken from density measurements. 

Other applications of the spectrometer’s automated features have been 
employed for other materials, particularly with superconductivity samples and 
organic materials. The superconductivity samples have included compounds con- 
taining Y, Ba, Sr, Cu, and 0. By modifying the compositions of the samples, 
the trends in infrared response were assessed along with a superconducting 
transition temperature that was unique to each compound. [14] 

The organic compounds have included succinonitrile (SCN) in glycerol; 
the spectrometer was used to assess the variations in infrared response in the 
different growth regions of this material . [15] 

The results above show that the spectrometer facility is capable of 
presenting data in numerous forms for analysis. However, the system is not 
restricted to the examples discussed. Spatial maps and other figures are 
available in color for easier interpretation, if desired, or compositions and 
other values may be plotted by number on a spatial map. Also, other optical 
parameters not mentioned previously (e.g., epilayer thickness derived from 
interference fringes) may be used for mapping. More variations for sample 
characterization will be innovated for the FTIR facility as new analytical re- 
quirements arise. 


CONCLUSION 

The FTIR spectrometer is a tool which has successfully demonstrated the 
capability of establishing the spatial property mapping of solid solution 
semiconductors. The results obtained from this system provide optical property 
characteristics as well as verify and complement the results from other 
characterization methods such as density measurements and x-ray dispersion 
analysis of the samples. In addition, the system provides a comprehensive, 
non-destructive analysis which is highly automated and takes a small amount of 
operator input and time. 
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Figure 2. Spatial composition map of a unidirectionally solidified sample of Hgj_ x Cd x Te cut perpendicular to 
the growth axis. The step size between points is 5 (250 ym), and each point represents the same area as the 
100 ym diameter aperture. 
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Figure 3. Spatial composition map of a unidirectionally solidified Hgj^CdjjSe wafer whose growth axis i 
horizontal from left to right. The step size between points is 10 (500 pm), and each point represents a 
300 ym diameter aperture size. 
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100 ym diameter aperture. 
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Figure 6. Contour map showing compositions of CdSe fraction, x, from a slice of 
Hgi- , x Cd x Se. The sample size and grid size are the same as in Figure 5. 
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HCT Grown w/Transverse Mag. Field 



13 


Figure 7. Plot showing the axial trend of CdTe composition, x, vs. distance of an ingot of Hg, Cd Te 
unidirectionally solidified using a transverse magnetic field. Both density measurement and FTlR data are 
shown. 




TABLE 1. FTIR HEADER FILE PARAMETERS 


1 apod triangle 

2 phasetype phaseapod 

3 xstart 4000.0 

4 xend 400.0 

5 scale none 

6 laseramp 2 

7 color red 

8 datasize w 

9 axescolor green 

10 datatype ras 

11 bkgscl 0 

12 sampscl 0 

13 speed 2 

14 ff symmetry double 

15 ymin 0.0 

16 ymax 1.0 

17 fvel 6 

18 rvel 6 

19 lpf 0 

20 hpf 0 

21 sampdir forward 

22 fstart 399.0 

23 fend 4001.0 

24 spec-range MIR 

25 dend 2548 

26 dstart 0 

27 nphzfft 512 

28 nphzdata 256 

29 nrawdata 2548 

30 nfft 4096 

31 ndata 2048 

32 nscans 256 

33 nscnbkg 16 

34 iris 50 

35 sgain 3 

36 title 

37 det 3 

38 material HgCdTe 

39 shape 0 

40 apert 100 

41 center 250 260 

42 radius 2.2 

43 transpl 1000.0 

44 transp2 1500.0 

45 trine 100.0 

46 opaquel 2000.0 

47 opaque2 2200.0 

48 transmin 0.0030 


apodization function (triangle or boxcar) 

phase correction (phase, real, imaginary or power) 

starting wavenumber for graph 

ending wavenumber for graph 

type of plot scaling (none or auto) 

sample every ith laser zero crossing (1 - 16) 

color for the plot 

size of data (w : 32 bit words s : 16 bit words) 
color for the axes 
default file extension 

scaling applied to .big file (not to be modified) 
scaling applied to .igm file (not to be modified) 
plotting speed (0-4 with 0 the fastest) 
type of fft symmetry (double or single) 
scaling for y-axis 

scaling for y-axis (yin = ymax for auto scaling) 

forward mirror velocity (1 - 7) 

reverse mirror velocity (1 - 7) 

low pass filter (0 - 8) 

high pass filter (0 - 2) 

sampling direction (forward, reverse, or both) 

starting wavenumber in file 

ending wavenumber in file 

spectrometer scanning range (FIR MIR NIR) 

ending point of an interferogram plot 

starting point of an interferogram plot 

number of phase fft points 

number of points in phase array 

number of input data points (nrawdata > ndata) 

number of fft points (must be power of 2) 

number of data points used 

number of scans to be co- added 

number of scans to be co- added for background 

iris opening size (%) 

post amplifier gain (1 - 3) 

’C2-H 2.40 cm’ 
detector selection (1 - 3) 
type of material 

0 = circular and 1 = rectangular 
aperture diameter in micrometers 
sample center in xy stage coordinates 
sample radius in cm 

wavenumber to begin searching for max. trans. 
wavenumber to stop searching for transmax 
wavenumber increment (dT) to average for transmax 
wavenumber to begin averaging for opaque transmittance 
wavenumber to stop averaging for opaque transmittance 
do not calculate cut-on if transmittancectransmin 


NOTE: Numbers 38-48 are the additional parameters used by the algorithm to extract 
cut-on wavelengths from transmission spectra. 
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